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Abstract 

oo 

qq This paper introduces QuanTree vl.l and QuanLin vl.l, two Java applications avail- 

able for free. (Source code included in the distribution.) Each application compiles a 
different type of input quantum evolution operator. The applications output a quan- 
^ turn circuit that is approximately equal to the input evolution operator. QuanTree 

compiles an input evolution operator whose Hamiltonian is proportional to the inci- 
> dence matrix of a balanced, binary tree graph. QuanLin compiles an input evolution 

operator whose Hamiltonian is proportional to the incidence matrix of a line (open 
string) graph. Both applications also output an error, defined as the distance in 
the Frobenius norm between the input evolution operator and the output quantum 
circuit. 



1 Introduction 

We say a unitary operator acting on a set of qubits has been compiled if it has been 
expressed as a SEO (sequence of elementary operations, like CNOTs and single-qubit 
operations). SEO's are often represented as quantum circuits. 
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There exist software (quantum compilers) like Qubiter[lJ for compiling arbi- 
trary unitary operators (operators that have no a priori known structure). This paper 
introduces two special purpose quantum compilers, QuanTree and QuanLin. They 
are special purpose in the sense that they can only compile unitary operators that 
have a very definite, special structure. 

QuanTree vl.O and QuanLin vl.l are two Java applications, available[2] for 
free. (Source code included in the distribution.) Each application compiles a different 
kind of input quantum evolution operator. The applications output a quantum circuit 
that is approximately equal to the input evolution operator. QuanTree compiles an 
input evolution operator whose Hamiltonian is proportional to the incidence matrix of 
a balanced, binary tree graph. QuanLin compiles an input evolution operator whose 
Hamiltonian is proportional to the incidence matrix of a line (open string) graph. 
Both applications also output an error, defined as the distance in the Frobenius norm 
between the input evolution operator and the output quantum circuit. 

Recently, Farhi-Goldstone-Gutmann (FGG) wrote a paper [3j that proposes a 
quantum algorithm for evaluating NAND formulas via a quantum walk over a tree 
graph connected to a line ("runway") graph. Their paper has inspired a flurry of 
papers expanding on their ideas. Among these papers is one [4] written by me, which 
provides all the theoretical underpinnings of QuanTree and QuanLin. Please refer 
to Ref.jl] and the source code of QuanTree and QuanLin if you have any technical 
questions that are no addressed here. To get a quantum circuit for the FGG algorithm 
requires first finding a quantum circuit for the evolution operators of a tree and 
line graph, which is what QuanTree and QuanLin do. A future paper will combine 
QuanTree and QuanLin software to give a quantum circuit for the full FGG algorithm. 

The standard definition of the evolution operator in Quantum Mechanics is 
U = e~ ltH , where t is time and H is a Hamiltonian. Throughout this paper, we will 
set t — — 1 so U = e lH . If H is proportional to a coupling constant g, reference to 
time can be restored easily by replacing the symbol g by —tg ) and the symbol H by 
-tH. 

2 QuanTree 

2.1 Input Evolution Operator 

A binary tree with A + 1 levels has 1 + 2 + 2 2 + . . . 2 A = 2 A+1 — 1 nodes (a.k.a. as 
states). To reach Ns = 2 A+1 states, we add an extra "dead" or "dud" node, labelled 
with the letter "d" . This d node is not connected to any other node in the graph. If 
we include this dud node, then the number of leaves Ni vs is exactly half the number 
of nodes: Ni vs = |2 A+1 = 2 A . We will often use N B for the number of bits and 
Ns = 2 Nb for the number of states. Therefore, N B = A + 1. For example, Fig|l] 
shows a binary-tree graph with A = 2 and N B = 3. It has Ns = 2 A+1 =8 nodes, 
labelled d ) 1, 2, . . . , 7, half of which (4, 5, 6, 7) are leaves. 
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4 5 6 7 

Figure 1: Binary tree with 8 nodes 



The Hamiltonian for transitions along the edges of the graph FigjTJis: 
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(1) 



where g is a real number that we will call the coupling constant. In Eq.Q, empty 
matrix entries represent zero. 

For N B = 3 qubits (i.e., 2 Nb = 8 states), the input evolution operator for 

QuanTree is U = e lHtree , where i^ ree is given by Eq. ([!]). It is easy to generalize Fig|l] 

AT_ M-.-.^TV™ ^^-^^^1^ AHtree £^ AT _ r fO Q /I 1 



and Eq.([lJ) to arbitrary 7V#. QuanTree can compile e lHtree for AT^ G {2, 3, 4, . . .}. 

For r — 1, 2, 3, . . if U = L r (g) + 0{g r+l ) ) we say L r (gf) approximates (or 
is an approximant) of order r for U . 

Given an approximant L r (g) + 0(g r+1 ) of £7, and some 7V T = 1, 2, 3 . . ., one 

/ \ Nt r+l 

can approximate !7 by f L r (^) J + O(f^y). We will refer to this as Trotter's trick, 
and to N T as the number of trots. 

For Nt = 1, QuanTree approximates e lHtree with an approximant of order 3 
that is derived in Ref.[4]. Thus, for N T = 1, the error is 0(g 4 ). For N T > 1, the error 
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2.2 The Control Panel 



Figj2] shows the Control Panel for QuanTree. This is the main and only window of 
the application. This window is open if and only if the application is running. 
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QuanTree 




Ver. 1.1 
Inputs— 

File Prefix 

Number of Qubifls 

Coupling Constant 

Number of Trots 

Order of Appro x. 



Output Highlights- 
Number of Elem. Ops 
Error 

Message 



i) 



(Write Files) 



Figure 2: Control Panel of QuanTree 



The Control Panel allows you to enter the following inputs: 

File Prefix: Prefix to the 3 output files that are written when you press the Write 
Files button. For example, if you insert test in this text field, the following 3 
files will be written: 

• test_qtree_log.txt (See Figj3]) 

• test_qtree_eng.txt (See Fig|3]) 

• test_qtree_pic.txt (See FigjS]) 

Number of Qubits: The parameter N B = 2, 3, 4, . . . defined above. 
Coupling Constant: The parameter jGR defined above. 
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Number of Trots: The parameter N T = 1, 2, 3, . . . defined above. 



The Control Panel displays the following outputs: 

Number of Elementary Operations: The number of elementary operations in 
the output quantum circuit. If there are no LOOPs, this is the number of lines 



in the English File (see Sec. 2.3.2), which equals the number of lines in the Pic- 
ture File (see Sec. |2.3.3| ). When there are LOOPs, the "LOOP k REPS : N T " and 
"NEXT k" lines are not counted, whereas the lines between "LOOP k REPS : N T " 
and "NEXT k" are counted N T times. 

Error: The distance in the Frobenius norm between the input evolution operator and 
the output quantum circuit (i.e., the SEO given in the English File). For a nice 
review of matrix norms, see Ref.[5j. For any matrix A G C nXn , its Frobenius 



norm is defined as \\A\\ F = y ^ k Aj^A* k . Another common matrix norm is 

the 2-norm. The 2-norm \\A\\2 of A equals the largest singular value of A. The 
Frobenius and 2-norm of A are related by[5]: \\A\\ 2 < \\A\\ F < V2\\A\\ 2 . Since 
the approximant used by QuanTree is of order 3, if e denotes the error, then 
e(g) « Kg 4 , for some K G R and \g\ « 1. Thus, 



log(e(g 2 )/e(g 1 )) _ ^ 
log(^M) 

For example, with Nb — 4 and Nt — 1, QuanTree gives e(0.05) = 1.383 x 10" 5 
and e(0.06) = 2.923 x 1(T 5 , which gives log(e 2 /ei)/log(^ 2 M) = 4.11. 

Message: A message appears in this text field if you press Write Files with a bad 
input. The message tries to explain the mistake in the input. 



2.3 Output Files 

Figsj3| [ij [5| were all generated in a single run of QuanTree (by pressing the Write 
Files button just once). They are examples of what we call the Log File, English 
File, and Picture File, respectively, of QuanTree. Next we explain the contents of 
each of these output files. 



2.3.1 Log File 

Figj3]is an example a Log File. The Log File records all the information found in the 
Control Panel. 
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© © © |3 test_qtree_log.txt 


[inputs : 


File Prefix = test 




Number of Bits = 3 




Coupling Constant =0.5 




Number of Trots = Z 




Ul Ucl UT ApfJf ox. Lrrlurl t Li u LWuyis J . 




Outputs : 




test_qtree_eng .txt 




test_qtree_pic .txt 




test_qtree_log .txt 




Number of Elem. Ops. = 15 




Error = S.lZ4S66e-03 









Figure 3: Log File generated by QuanTree 



2.3.2 English File 

Figj4]is an example of an English File. The English File completely specifies the out- 
put SEO. It does so "in English", thus its name. Each line represents one elementary 
operation, and time increases as we move downwards. 

In general, an English File obeys the following rules: 



• Time grows as we move down the file. 

• Each row corresponds to one elementary operation. Each row starts with 4 
letters that indicate the type of elementary operation. 

• For a one-bit operation acting on a "target bit" a ) the target bit a is given after 
the word AT. 

• If the one-bit operation is controlled, then the controls are indicated after the 
word IF. T and F stand for true and false, respectively. aT stands for a control 
n(a) at bit a. aF stands for a control n(a) at bit a. 

• "LOOP k REPS : 7V T " and "NEXT k" mark the beginning and end of N T Trotter 
iterations, k labels the loop, k also equals the line-count number (first line is 
0) of the line "LOOP k REPS:iV r " in the English file. 

• SWAP a (3 stands for the swap (exchange) operator E{a ) (3) that swaps bits a 
and (3. 
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eoe 



■ test_qtree_eng.txt 



|WTY 


-90.0 


AT 


IF 


ZT 








SWAP 


1 


IF 


ZT 












ROTY 


-90.0 


AT 


I r 


IT 










LOOP 


3 


REPS : Z 














ROTX 


ZO.Z57117113534SS6 


AT 


1 


I 7 


Z r 






ROTX 


40.3O355Z05O94Z35 


AT 


Z 


IF 


IT 


0T 




ROTX 


-Z0. 257117113 5 348 8 6 


AT 
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IF 


zr 






ROTX 


40.936Z575OOZ6S4Z 


AT 


L 


T — 

IF 


z- 


OT 




ROTX 


4O.514Z34ZZ706978 


AT 


Z 


IF 


IT 


OF 




NEXT 


3 
















ROTY 


90.0 


AT 


I r 


IT 










SWAP 


1 


IF 


ZT 












ROTY 


90.0 


AT 


I 7 


ZT 











Figure 4: English File generated by QuanTree 



• PHAS stands for a controlled one-bit gate, where the one-bit gate consists of 
diag{\^ 1) times an angle ("phase"). 

• POPH stands for a controlled one-bit gate, where the one-bit gate consists of 
Po = n() times an angle ("phase"). P1PH stands for a controlled one-bit gate, 
where the one-bit gate consists of Pi = n() times an angle ("phase"). 

• SIGX, SIGY, SIGZ, HAD2 stand for the Pauli matrices dj, oy, az and the one-bit 
Hadamard matrix H. 

• ROTX, ROTY, ROTZ, ROTN stand for rotations with rotation axes in the directions: 
x, z, and an arbitrary direction n, respectively. 

Here is a list of examples showing how to translate the mathematical notation 
used in Ref.[4j into the English File language: 
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eee ii 


test_qtree_pic.txt 


^— +— Ry 






P— <— > 






1 e>— Ry 






LOOP 3 


REPS : Z 




0— Rx 1 






Rx— P— P 






0— Rx 1 






0— Rx— & 






Rx— P— 






NEXT 3 






1 fi>— Ry 






§>—<—> 






P— +— Ry 




/a 



Figure 5: Picture File generated by QuanTree 



Mathematical language 


English File language 


Loop called 5 with 2 repetitions 


LOOP 5 REPS: 2 


Next iteration of loop called 5 


NEXT 5 


£( 1;0 )n(3W2) 


SWAP 1 IF 3F 2T 


e i42.7n(3)n(2) 


PHAS 42.7 IF 3F 2T 


e i42.7n(3)n(2) 


POPH 42.7 AT 3 IF 2T 


e i42.7n(3)n(2) 


P1PH 42.7 AT 3 IF 2T 


(7jc (l)n(3)n(2) 


SIGX AT 1 IF 3F 2T 


ay (l)n(3)n(2) 


SIGY AT 1 IF 3F 2T 


az (l)n(3)n(2) 


SIGZ AT 1 IF 3F 2T 


#(l)™(3)n(2) 


HAD2 AT 1 IF 3F 2T 


( e |if 5 23.7 ( T X (l)^n(3)n(2) 


ROTX 23.7 AT 1 IF 3F 2T 


( e |ifo23.7 ( r y (l)^n(3)n(2) 


ROTY 23.7 AT 1 IF 3F 2T 


( e ilf(5 23 - 7<T 2(l))n(3)n(2) 


ROTZ 23.7 AT 1 IF 3F 2T 


( e ilf(5[ 30<T ^( 1 )+ 40<T y( 1 )+ 11 ^(l)])n(3)n(2) 


ROTN 30.0 40.0 11.0 AT 1 IF 3F 2T 
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2.3.3 ASCII Picture File 



Figj5] is an example of a Picture File. The Picture File partially specifies the output 
SEO. It gives an ASCII picture of the quantum circuit. Each line represents one 
elementary operation, and time increases as we move downwards. There is a one-to- 
one onto correspondence between the rows of the English and Picture Files. 
In general, a Picture File obeys the following rules: 



• Time grows as we move down the file. 

• Each row corresponds to one elementary operation. Columns 1, 5, 9, 13, . . . rep- 
resent qubits (or, qubit positions). We define the rightmost qubit as 0. The 
qubit immediately to the left of the rightmost qubit is 1, etc. For a one-bit 
operator acting on a "target bit" a, one places a symbol of the operator at bit 
position a. 

• I represents a wire connecting the same qubit at two times. 

• -represents a wire connecting different qubits at the same time. 

• + represents both I and -. 

• If the one-bit operation is controlled, then the controls are indicated as follows. 
@ at bit position a stands for a control n(a). at bit position a stands for a 
control n(a). 

• "LOOP k REPS : A T " and "NEXT k" mark the beginning and end of N T Trotter 
iterations, k labels the loop, k also equals the line-count number (first line is 
0) of the line "LOOP k REPS: A T " in the picture file. 

• The swap (exchange) operator E(a, (3) is represented by putting arrow heads < 
and > at bit positions a and (3. 

• A phase factor e l ° for some angle 9 is represented by placing Ph at any bit 
position which does not already hold a control. 

• The one-bit gate Po(ct) = n(a) times an angle is represented by putting OP at 
bit position a. 

• The one-bit gate Pi (a) = n(a) times an angle is represented by putting @P at 
bit position a. 

• One-bit operations o"x(<^), cry (a), (Jz{ol) and H(a) are represented by placing 
the letters X,Y,Z, H, respectively, at bit position a. 
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• One-bit rotations acting on bit a, in the x,y,z,n directions, are represented by 
placing Rx,Ry,Rz, R, respectively, at bit position a. 

Here is a list of examples showing how to translate the mathematical notation 
used in Ref.[4j into the Picture File language: 



Mathematical language 


Picture File language 


Loop called 5 with 2 repetitions 


LOOP 5 REPS: 2 


Next iteration of loop called 5 


NEXT 5 


V i ) 


@ < > 


e i42.7n(3)n(2) 


@ + — Ph 


e i42.7n(3)n(2) 


OP — @ I I 


e i42.7n(3)n(2) 


@P — © I I 




@ X | 


(1) n(3)n(2) 


@ Y | 


az (l)n(3)n(2) 


0— @— Z | 


^Q\n(3)n(2) 


@ H | 


( e j 1 fo23.7a x (l)jn(3)n(2) 


@ Rx | 


^§^23.7^(1)^(3)71(2) 


@ Ry | 


^§^23.7^(1)^(3)72(2) 


@ Rz | 




© R 1 



2.4 Behind the Scenes 

Brief summary of the steps taken by QuanTree every time you press the Write Files 
button: 

1. Generate the English and Picture Files according to the rules of Ref.[4j. 

2. Generate H . Calculate the eigenvalues and eigenvectors of H. Use this eigen- 
system to calculate U = e lH . 

3. Read the English File that was written in Step[T} Multiply out the SEO given by 
the English File to obtain a unitary matrix U f . Calculate the error \\U — U f \\p. 

4. Generate the Log File. 
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3 QuanLin 

3.1 Input Evolution Operator 




000 001 010 Oil 100 101 110 111 




Figure 6: Line (open string) with 8 nodes 



Let Nb be the number of bits and N$ = 2 Nb the number of states. Fig |6] shows 
the 8 possible states for 3 bits. The states read from left to right are in increasing 
"decimal ordering" . These 8 states can also be ordered in "Gray code ordering" . A 
sequence of Gray code is one wherein two consecutive states are labelled by a binary 
number and these labels differ only at one bit position. In FigJHJ states connected by 
an edge (curved line) are consecutive in a Gray code ordering. 

The Hamiltonian for transitions along the edges of the graph FigJHjis: 
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(3) 



where g is a real number that we will call the coupling constant. In Eq.Q, empty 
matrix entries represent zero. 

For Nb = 3 qubits (i.e., 2 Nb = 8 states), the input evolution operator for 
QuanLin is U — e lHline , where Hu ne is given by Eq.Q. It is easy to generalize Figj6] 
and Eq.Q to arbitrary Nb- QuanLin can compile e lHline for Nb g{2,3,4,...}. 

For r = 1, 2, 3, . . ., if U = L r (g) + 0{g r+l ) ) we say L r (g) approximates (or 
is an approximant) of order r for U . 



11 



Given an approximant L r (g) + 0(g rJrl ) of [7, and some N T = 1, 2, 3 . . one 

can approximate J7 by ( L r (^) J + C?(^f). We will refer to this as Trotter's trick, 

and to N T as the number of trots. 

For Nt = 1, QuanLin approximates e lHline with a Suzuki approximant of order 
r = 2, 4, 6, . . . that is derived in Ref.[4j. Thus, for Nt = 1, the error is 0(g r+1 ). For 
7V T > 1, the error is O(^). 

3.2 The Control Panel 

Fig{7] shows the Control Panel for QuanLin. This is the main and only window of 
the application. This window is open if and only if the application is running. 



eoe 



QuanLin 




Ver. 1.1 
Inputs— 

File Prefix 

Number of Qubits 

Coupling Constant 

Number of Trots 

Order of Approx. 



T 



ID 



J3 



' T J 



Output Highlights- 
Number of Elem, Ops 
Error 

Message 



(Write Files) 



Figure 7: Control Panel of QuanLin 



The Control Panel for QuanLin is almost identical to that for QuanTree. The 
significance of the various data fields in the Control Panel for QuanLin is the same 
as for QuanTree. 
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3.3 Output Files 



Figsj8j [9j [lOj were all generated in a single run of QuanLin (by pressing the Write 
Files button just once). They are examples of what we call the Log File, English 
File, and Picture File, respectively, of QuanLin. These files are analogous to their 
namesakes for QuanTree. They follow the same rules. 



8 G e\ test_qline_log.txt 


[inputs : 


File Prefix = test 




Number of Bits = 3 




Coupling Constant =0.5 




Number of Trots = Z 




Order of Suzuki Approximant = 4 




Outputs : 




test_qline_eng . . tstt 




test_qline_pic .tsct 




test_qline_log .tsct 




Number of El em. Ops. = 3Z 




Error = Z.9Z6344e-04 









Figure 8: Log File generated by QuanLin 
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©ee 


■ ■ toct nil 


i r ic_c r i y * i a. l 






LOOP 


REPS: Z 










KUI A 


D . =r 3 f 14£ =rb /Y 344bH 


Al 


/■> 

<v 








KUI a 


11 . B r 4£HbyjD4bB3J r 


AT 


1 IF 


0T 






ROTX 


11.874285935468937 


AT 


Z IF 


IT 


0F 




ROTX 


11.S74ZS5935468937 


AT 











ROTX 


11.S74ZS5935468937 


AT 


1 IF 


0T 






ROTX 


11.374Z35935463937 


AT 


Z IF 


IT 


0F 




ROTX 


-3.4874840Z493Z8ZZ7 


AT 










ROTX 


-1S.849Z5398533458 


AT 


1 IF 


0T 






ROTX 


-18 . 849Z5398533458 


AT 


Z IF 


IT 


0F 




ROTX 


-3.4874S40Z493ZSZZ7 


AT 











ROTX 


11.874285935468937 


AT 


1 IF 


0T 






ROTX 


11.S74ZS5935468937 


AT 


Z IF 


IT 


0F 




ROTX 


11.874Z85935468937 


AT 











ROTX 


11.S74Z3593546S937 


AT 


1 IF 


0T 






ROTX 


11.S74Z85935468937 


AT 


Z IF 


IT 


0F 




ROTX 


5.93714Z967734468 


AT 











NEXT 






























Figure 9: English File generated by QuanLin 
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O f> 

Vj* 'S—r' V_r* 


i) test_qline_pic.txt 




REPS: Z 




1 1 Rx 

1 I It. A. 






1 Rx— 






Rx— ^ — 






1 1 Rx 






1 Rx— 






Rsc— §— 






1 1 Rx 






1 Rx— & 






Rx— 0— 






1 1 Rx 






1 Rx— 






Rx— P— 






1 1 Rx 






1 Rx— 






Rx— 






1 1 Rx 






NEXT 




/a 



Figure 10: Picture File generated by QuanLin 
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